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Abstract. The multiplicity (iVch) and pseudorapidity distribution {dNch/drj) of 
primary charged particles in p+p collisions at Large Hadron Collider (LHC) energies 
of y/s — 10 and 14 TeV are obtained from extrapolation of existing measurements 
at lower -y/s. These distributions are then compared to calculations from PYTHIA 
and PHOJET models. The existing ^/s measurements are unable to distinguish 
between a logarithmic and power law dependence of the average charged particle 
multiplicity ((TVch)) on ^/s, and their extrapolation to energies accessible at LHC give 
very different values. Assuming a reasonably good description of inclusive charged 
particle multiplicity distributions by Negative Binomial Distributions (NBD) at lower 
^/s to hold for LHC energies, we observe that the logarithmic y/s dependence of 
(Nch) are favored by the models at midrapidity. The dNch/drj versus rj for the 
existing measurements are found to be reasonably well described by a function with 
three parameters which accounts for the basic features of the distribution, height at 
midrapidity, central rapidity plateau and the higher rapidity fall-off. Extrapolation 
of these parameters as a function of ^/s is used to predict the pseudorapidity 
distributions of charged particles at LHC energies. dNch/dr] calculations from PYTHIA 
and PHOJET models are found to be lower compared to those obtained from the 
extrapolated dNch/drj versus rj distributions for a broad rj range. 
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1. Introduction 

The Large Hadron Collider (LHC) at CERN is designed for colliding proton-proton 
beams upto ^/s = 14 TeV [Ij. Collisions at these unprecedented high energies will 
provide opportunities for new physics [Ij. In order to fully exploit the enormous 
physics potential it is important to have a complete understanding of the reaction 
mechanism. The particle multiplicity distributions, one of the first measurements to 
be made at LHC, will be used to test various particle production models based on 
different physics mechanism and also provide constrains on model features. Some of 
these models are based on string fragmentation mechanism p] and some are based on 
Pomeron exchange [3j. 

In this paper, we first make a compilation of the existing data on the average 
charged particle multiplicity (at midrapidity and full rapidity range) and charged 
particle pseudorapidity distribution as a function of ^/s. Then we judiciously extrapolate 
the measurements to obtain prediction of {Nch) (for midrapidity and full rapidity range), 
charged particle multiplicity distributions and dNch/drj distributions at LHC energies 
of a/s = 10 and 14 TeV. These results are also compared to calculations from PYTHIA 
\2\ and PHOJET [3] monte carlo models. Energy dependence of average number of 
particles produced in p+p{p) collisions can be used to distinguish various models [1]. 
A statistical model [5] and some hydrodynamical models |6j predict a dependence of 
{N) ~ s^/^, whereas multiperipheral models [7j and Feynman's scaling [8] lead to the 
dependence as (A^) ~ ln{s). Such a logarithmic dependence and on higher powers of 
ln{s) is also predicted by Regge- Mueller model [9]. Arguments based on simple phase 
space considerations [10] however predict a power law dependence as (A^) ~ s^^'^. In this 
work we show that for the (A"ch) data available at various ^/s, we cannot distinguish 
the logarithmic and power law dependences on ^/s. However the measurements at 
the LHC energies will provide a clear answer. The pseudorapidity distributions on 
the other hand is found to be described by a form which resembles a generalized 
Fermi distribution. Such distributions have been used to explain the pseudorapidity 
distributions of produced particles in hadronic collisions at ISR [11] and heavy-ion 
collisions at RHIC [12]. 

The compiled experimental data presented in this paper corresponds to Non Singly 
Diffractive (NSD) events for minimum bias triggers. The charged particle data were 
corrected for secondary interactions, gamma conversions, short lived decays {Kg, A), 
reconstruction efficiency and acceptance effects by the experiments. To match the 
experimental conditions, the model simulations presented are also corrected for short 
lived decays. A transverse momentum [px) cut of greater than 100 MeV/c are usually 
used in realistic experimental conditions, as for example in ALICE experiment at LHC. 
Model simulations using PYTHIA and PHOJET suggest, a 6% ± 2% effect on the 
charged hadron multiplicity due to a 100 MeV/c cut-off in px at y/s = 10-14 TeV. The 
error of 2% comes from the difference in results from PYTHIA and PHOJET models. 
At LHC, while most of the experiments will have mid-rapidity measurements of charged 
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particle multiplicity, ALICE experiment has the possibility to measure the distributions 
over -5.0 < rj < 3.5 range [1]. The CMS and ATLAS experiments will have a more 
limited coverage of |?7| < 2.5 units [1]. 

The model results presented are from PYTHIA using version 6.4 (ATLAS tuned) 
and those from PHOJET with version 1.12 (default settings). It may be mentioned 
that recently a new tuned version of PYTHIA has been released [13]. The PHOJET 
model combines the ideas based on a dual parton model [H] on soft process of 
particle production and uses lowest-order perturbative QCD for hard process. Regge 
phenomenology is used to parameterize the total, elastic and inelastic cross-sections. 
The initial and final state parton shower are generated in leading log- approximation. 
PYTHIA on the other hand uses string fragmentation as a process of hadronization 
and tends to use the perturbative parton-parton scattering for low to high px particle 
production. Although there are several other theoretical predictions on total cross 
section expected at LHC energies [15], current work focuses on how a judicious 
extrapolation from existing multiplicity data compares to the calculations from some of 
the available models [21 [3]. 

2. Multiplicity distribution 




Figure 1. Multiplicity distribution for charged particles in p+p collisions at various 
center of mass energies at midrapidity (left panel) and full rapidity (right panel) 
ranges [16]. The errors are statistical. The solid lines are NBD fit to the data points 
using the function given in Eqn. [1] 



The measurements of charged particle multiplicity distribution has been found to 
be well described by negative binomial distribution (NBD) at midrapidity and also for 
the full rapidity region in p+p p!6] . The NBD distribution has a form. 
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Figure 2. NBD fit parameters (A'ch) (left panel) and 1/fc (rigfit panel) to multiplicity 
distributions in p+p collisions at various center of mass energies |16| . The lines are 
fits to the data points using the function forms as discussed in text. Also shown for 
comparison are the results from PYTHIA and PHOJET models at ^/s = 200, 540 and 
900 GeV. 



The NBD has two parameters, (A^^ch) and k. Where the parameter k is an interesting 
quantity, would correspond to Poisson distribution (independent particle 

production) and k = 1 would correspond to Geometric distribution. Under the limit 
of large multiplicity {Net — ^ Large), the NBD distribution goes over to a Gamma 
Distribution. Some of the measured multiplicity distributions at midrapidity (|?7| < 0.5) 
and over full pseudorapidity range (|?7| < 5.0) are fitted to NBD distribution and are 
shown in Fig. [1] left and right panels respectively. The NBD parameters (A^ch) and k 
extracted are plotted in Fig. [2] left and right panels respectively. The CDF experiment 
results are not used in the current studies as they have multiplicity distributions with 
a high pt cut-off of 400 MeV/c measured within \r]\ < 1 [17J. PYTHIA and PHOJET 
model calculations suggest that there is more than 50% loss in the average number of 
charged particles due to a pr cut-off of 400 MeV/c at midrapidity. Since all other results 
presented have a much smaller pt cut-off ~ 100 MeV/c, inclusion of CDF results would 
make this comparative study heavily dependent on the model based extrapolation to 
lower Pt regions. 

It is observed from Fig. [2] that both at midrapidity and full rapidity range the 
(Nch) increases with ^/s while the k value decreases with ^/s. The ^/s region for which 
measurements exists, the (Nch) dependence can be described reasonably well by both 
the following expressions (as shown in Fig. [2]), 

(Ar,h) = a + b ln{^s) + c {ln{^s)f (2) 

and 

(iVch) = a' + b' {^sr\ (3) 

where a, 6, c, a', b' and d are fit parameters. The values of a, 6, c, a', b' and d at 
midrapidity are 3.8 ± 0.1, 1.5 ± 0.2, 0.45 ± 0.1, 2.3 ± 0.14, 1.5 ± 0.13 and 1.2 ± 0.26 
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respectively. For full rapidity coverage the values of a, b, c, a', U and c' are 1.3 ± 0.3, 
0.62 ± 0.17, 0.59 ± 0.02, -10.5 ± 0.86, 9.9 ± 0.69 and 0.22 ± 0.01 respectively. It is 
noted that extrapolation of the power law function to LHC energies seems to lead to 
a sudden increase in average charged particle multiplicity at midrapidity. Such a large 
unexpected jump in the multiplicities at midrapidity already seems to put constraints on 
applicability of such a functional form. As will be seen later in Fig. H] such a functional 
form will lead to a probability distribution of charged particle multiplicity at midrapidity 
not showing the characteristic drop at higher multiplicity. Such a drop is expected 
from PYTHIA, PHOJET models and is seen for available experimental data shown in 
Fig.lll In Fig. [2]the average number of charged particles from experiments at ^/s = 200, 
540 and 900 GeV are compared to corresponding results from PYTHIA and PHOJET 
models. Both the models seem to be in reasonable agreement with the measurements 
at midrapidity and full rapidity. 

The 1/k dependence on ^/s has the form a + f3 ln{^/s), where a and f3 are fit 
parameters. For midrapidity the values of a and (3 are 0.65 ± 0.03 and 0.07 ± 0.03 
respectively. For full rapidity case, the values of a and (3 are -0.11 ± 0.006 and 0.06 ± 
0.001 respectively. The results from PYTHIA and PHOJET models are also shown. At 
higher energies the results from PHOJET is in better agreement with the measurements 
compared to those from PYTHIA. For the midrapidity measurements PHOJET model 
calculations fail to match the data for -^i = 200 GeV. 

We have examined the charged particle multiplicity distribution from E735 
Collaboration for p+p collisions at y/s = 1800 GeV [18] . The multiplicity distribution 
is found to be well explained using the sum of two NBD functions instead of a single 
NBD as for the rest of the data discussed above. The sum of two NBD is given as, 

F = ujP^^MNci.)i, h; n) + (1 - cu)P^BD((A^ch)2, ^2; n) (4) 

where -P^bd -^nbd have the same form as in Eqn [H u; is the weight factor, 
(A'ch)i and {Nch)2, ^1 and k2 are the respective NBD parameters. The fit to the E735 
data by this function is shown in Fig. [31 The values of various parameters obtained from 
the fit are, (iVch)i = 36.5 it 1.7, (A/'ch)2 = 86.9 ± 2.6, h = 2.8 ± 0.2, ^2 = 10.3 ± 0.8 
and uj = 0.18 ± 0.03. The overall average value of the charged particle multiplicity is 
44.4. Since this data could not be fitted to a single NBD, in contrast to the data from 
other energies discussed in this paper, we have not included this measurement in our 
extrapolation studies. Further it is observed from the Fig. [2] that the two extrapolations 
discussed in Eqns [2] and [3] tend to differ only above y/s = 2000 GeV. The physics reason 
attributed to a different shape of multiplicity distribution at -^i = 1800 GeV compared 
to those at lower energies is due to multiple parton interactions |18]. These result are 
also indicative of the deviation from the KNO (Koba, Nielsen and Olesen) scaling |19j . 
The weight factor shows the second NBD distribution dominates among the two. It 
will certainly be interesting to find out if this feature is more pronounced at the LHC 
energies. 

The extrapolated values of {Nch) and k at both midrapidity and full rapidity regions 
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Figure 3. Relative cross-section of charged particles produced in p+p collisions at ^/s 
= 1800 GeV [18] fitted to a sum of to NBD functions. See text for more details. 



Table 1. Extrapolated NBD parameters (A^ch) and k for different y/s at midrapidity 
and full rapidity range for p+p collisions at a/s — 10 and 14 TeV. 



Eqn. 


(TeV) 


(iVch) 


k 


f] range 


Eqn. 


m 


10 


9.6 ± 1.1 


1.2 ± 0.14 


1^1 


< 0.5 


Eqn. 


m 


10 


56.7 ± 3.8 


2.2 ± 0.09 


1^1 


< 5.0 


Eqn. 


m 


10 


25.0 ± 11.5 


1.24 ± 0.14 


1^1 


< 0.5 


Eqn. 


m 


10 


63.7 ± 11.3 


2.2 ± 0.09 


1^1 


< 5.0 


Eqn. 


m 


14 


10.8 ± 1.3 


1.20 ± 0.14 


1^1 


< 0.5 


Eqn. 


m 


14 


60.7 ± 4.0 


2.14 ± 0.09 


1^1 


< 5.0 


Eqn. 


m 


14 


36.2 ± 18.5 


1.20 ± 0.14 


1^1 


< 0.5 


Eqn. 


m 


14 


69.4 ± 12.2 


2.14 ± 0.09 


1^1 


< 5.0 



using the functional forms as in Eqns [2] and [3] for ^/s = 10 and 14 TeV are given in 
the Table [H Note the difference in values of (Nch) depending on the dependence on a/s 
as per Eqn. [2] (logarithmic dependence) or Eqn. [3] (power law dependence). Knowing 
these values (parameters of NBD function) we can predict the multiplicity distributions 
for both midrapidity and full rapidity ranges at y/s = 10 and 14 TeV using Eqn. [TJ 
These distributions for both midrapidity regions and full rapidity regions are shown in 
Fig. m and Fig. [5] respectively. The left panels are for y/s = 10 TeV and right panels 
are for ^/s = 14 TeV. The results when compared to PYTHIA and PHOJET model 
calculations show that the extrapolation of (A^^ch) using Eqn. [2] is favored by the models 
at midrapidity. The extrapolated results show no such preference to models in full 
rapidity region. In general PYTHIA results are found to be higher compared to those 
from PHOJET calculations. Actual experimental measurements at LHC will confirm 
the preferred (Ach) dependence on y/s. 
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Figure 4. Estimated multiplicity distribution for charged particles in p+p collisions 
at y/s = 10 (left panel) and 14 TeV (right panel) in midrapidity. Solid and dashed lines 
are distributions obtained from (A^ch) extrapolation using Eqns.[2]and [3] respectively. 
The dotted lines reflects errors in multiplicity distributions due to extrapolation of the 
parameters (Nch) and k. The results are compared to corresponding calculations from 
PYTHIA and PHOJET. 




600 



Figure 5. Same as Fig. |4] for full rapidity region. 
3. Pseudorapidity distribution 

The pseudorapidity distributions of charged particles from the existing data at ^/s = 
53, 200, 546 and 900 GeV [20] can be described by the following functional form, 

dN _ C + r] 

Id^ ~ 1 + exp» ^ ' 

This formula is chosen to describe the central plateau and the fall off in the 
fragmentation region of the distribution by means of the parameters ?7o and 5 
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Figure 6. Left panel: Pseudorapidity distribution for charged particles in p+p 
collisions at various center of mass energies . The solid lines are fit to the data points 
using the function given in Eqn. [5l Right panel: Comparison of the pseudorapidity 
distributions for charged particles at various ^/s to PYTHIA and PHOJET model 
calculations. 



Table 2. Parameters C, rjo and 6 for different -^/s. 



Collision 


(GeV) 


C 


^0 


6 


p+p 


53 


2.4 ± 0.23 


1.5 ± 0.17 


1.12 ± 0.1 


p+p 


200 


2.5 ± 0.07 


2.5 ± 0.05 


1.10 ± 0.04 


p+p 


546 


3.0 ± 0.10 


2.9 ± 0.07 


1.15 ± 0.04 


p+p 


900 


3.6 ± 0.10 


3.0 ± 0.05 


1.36 ± 0.05 



respectively. The term C + rj, C is a fit parameter, describes the magnitude of the 
distribution and the dip at the 77 = 0. A similar form has been used to describe the p+p 
data at ISR energies [H] and heavy- ion collisions at RHIC [12]. The distribution is a 
generalization of Fermi distribution and recent work suggests a relation of this functional 
form to string percolation model |21j . 

The values of the parameters C, r/o and 6 obtained by fitting the data distributions 
with Eqn. [5] are given in Table [2] and the fits to data are shown in Fig. [6] (left 
panel). The value of parameters C and rjQ are found to increase with increasing 
^/s. The value of the parameter 6 is found to be approximately independent of a/s 
within errors. The constancy of 6 is another way of demonstrating the concept of 
limiting fragmentation [22]. In such a scenario, multiplicity density in pseudorapidity 
when plotted as a function of pseudorapidity shifted by beam rapidity is expected 
to be independent of pseudorapidity at forward rapidities p2j. Shown in the right 
panel of Fig. [6] are the comparisons of the experimentally measured charged particle 
pseudorapidity distributions to those from the PYTHIA and the PHOJET calculations. 
The CDF measurements [23] are not used to obtain the parameters C, tjq and S for 
making predictions at LHC energies as discussed later, because of their very limited r) 
coverage. It is observed that for the lowest beam energy studied, -y/i = 53 GeV, both 
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Figure 7. Expected pseudorapidity distribution (solid lines) for charged particles in 
p+p collisions at y/s = 10 (left panel) and 14 TeV (right panel). This is obtained 
from the extrapolation using the existing data from p+p collisions at lower energies. 
The dotted lines indicate the uncertainties associated with the extrapolation. Also 
shown for comparison are the expected dNch/drj from PYTHIA and PHOJET model 
calculations for p+p collisions at \/s = 10 and 14 TeV. 



PYTHIA and PHOJET are in very good agreement with the experimental data. As the 
beam energy increases PYTHIA results seems to be in better agreement with the data. 
While for the top energy {y/s = 1800 GeV |17j) studied the error bars are large to make 
a conclusion on which of the two models have better agreement with the data. 

Using the average value of 6 and extrapolating the value of C and r^o to y/s = 10 and 
14 TeV, we are able to predict the full pseudorapidity distribution for charged particles. 
The extrapolation is done by fitting the variation of C with ^/s with a functional form 



The values of C obtained are 7.63 ± 0.87 and 8.43 ± 1.04 for = 10 and 14 TeV 
respectively. Those for rjo are 4.01 ± 0.14 and 4.15 ± 0.15 for ^/s = 10 and 14 TeV 
respectively. Assuming the functional form given in Eqn. [5] is valid for ^/s = 10 and 
14 TeV and using the parameter values obtained as above, we can now predict the 
full pseudorapidity distribution of charged particles at ^/s = 10 and 14 TeV. These 
distributions are shown in Fig. [7] along with expectations calculated from PYTHIA and 
PHOJET models for the same collisions. The left panel shows the results for ^/s = 10 
TeV and right panel shows the results for ^/s = 14 TeV. 

It is observed that at both y/s = 10 and 14 TeV, the predictions from a judicious 
extrapolation of existing data in general are above the model predictions. The PYTHIA 
results are close to the lower error band (dotted lines) for r/ < 4. The values 
from PYTHIA are higher than those from PHOJET model calculation over a large 



C = 3.7 + 1.15 ln{^/s) + 0.25 iln{^/^)Y 

and the variation of tiq with y/s as 

rio = 3.1 + 0.4 /n(v^). 



(6) 



(7) 
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pseudorapidity range presented. In general all the three distributions have almost similar 
shape. The differences between the models could arise due to several reasons. The event 
generation in PYTHIA is mainly designed to describe the possible hard interactions in 
p+p{p) collisions. It also combines sophisticated models dealing with soft hadronic 
interactions [2]. However in the case of PHOJET, the main approach is to describe 
the soft component of hadron-hadron, photon-hadron interactions at high energies. It 
combines the hard component calculated by perturbative QCD at the partonic level [3]. 
Due to the different underlying theoretical models used PYTHIA has larger model 
parameters which can be adjusted in order to better reproduce the data compared 
to PHOJET. Another difference lies in the parametrisation used get the p+p{p) cross- 
sections. PYTHIA uses those derived from the Pomeron exchange model while the 
PHOJET uses the optical theorem and cross-sections are corrected for high energies 
using the unitarity principle. 

4. Summary 

We have obtained the (Nch), distribution and dNch/drj versus rj for ^/s = 10 and 
14 TeV using extrapolations of existing data. The results have been compared to cal- 
culations from PYTHIA and PHOJET. Measurements at midrapidity at LHC will help 
distinguish whether (Nch) has a power law or a logarithmic dependence on -y/s. The 
shape of dNch/dr] versus t] obtained from the extrapolation of information from existing 
measurements are very similar to those calculated from models. PYTHIA calculations 
of A'ch and dNch/drj distributions seem to be higher compared to those from PHOJET 
model. The N^h distributions from models suggest the multiplicity distribution at LHC 
energies may not be well described by a single NBD distribution, already such a be- 
haviour has been seen at ^/s = 1.8 TeV [T8] . 
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